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HIERARCHICAL MARKOVIAN ALGORITHM IN QCD EVOLUTION " 



A. KUSINA 

Institute of Nuclear Physics, Polish Academy of Sciences, 

ul. Radzikowskiego 152, 31-34-2 Cracow, Poland. 

A general formalism was introduced for reorganization of QCD evolution equations and deriva- 
tion of hierarchical solution to DGLAP equation. The hierarchical solution separates two types 
of parton emissions: the flavour changing emissions and the diagonal ones (bremsstrahlung). 
It has been shown that both of the achieved processes are of the Markovian type, what gives a 
nice possibility of Monte Carlo implementation. Hierarchical algorithm has been implemented 
and crosschecked with standard Markovian algorithm implemented in program EvolFMC. 



1 Introduction 



Qh General Form of QCD evolution equation is given by 

& dtDk{t,x)=Y,Pkj{t,-)®D,{t,-){x), (1) 

, where the convolution rule is defined by 



^ /l(-) » /2(-)(^) = j dZidZ26{x - Z^Z2)fl{zi)f2{z2). (2) 

CO Parton momentum distribution function (PDF) Dk is a function of fraction of hadron longitu- 

■^ dinal momentum - x and evolution variable t = InQ (where Q is virtuality). Splitting function 

^^ Pkj which represents probability of parton emission, can be calculated perturbatively. General 



form of Pi-i is 



> 



:>& l^-\fi..x(^ ^\ , zaO 



Pu,{t, Z) = -P^^{t)6,k6{l -z) + P^^it, Z). (3) 



KA We impose the momentum sum rule: X^fc /o dx xDk{t, x) = const what leads to relation between 

PLit)=J2fdzzP^,{t,z). (4) 



5—1 virtual and real part of kernel 



J '" 



2 Reorganisation of evolution equation 

Our aim is to separate two types of parton emissions: flavour changing emissions and the diagonal 
ones (bremsstrahlung) . To do it we divide kernel into two parts the diagonal in flavour part and 
the rest 

Pf,{t,z) = -P^,^{t)6,k5{l -z) + 6,kPg{t,z), 

(5) 
Pf.it, z) = P.kit, z) - Pf,{t, z) = -Pf,^(t) V<5(1 - ^) + (1 - 5,k)P%{t. z). 

■^Talk given at XLIIId - Rencontres de Moriond QCD and High Energy Interactions, La Thuile, Italy, March 
08-15 2008. 



We do it in such a way that for both parts we have relation analogous to equation Q. We 
will use the compact matrix notation, where k and x are treated as indices. But we need to 
remember that we do not have the properties of commutativity anymore. Now the evolution 
equation has the form: 



dMt) = P(t)D(t) = (P^(t) + P^(t))D(t). 



(6) 



We introduce the diagonal transitions operator Gyi(t, io) which is solution of the equation ana- 
logues to our evolution equation ([I]) but with only diagonal part of the kernel: 



dtGA{t,to)=P'^{t)GA{t,to). 
Solution to this equation can be written as time ordered exponent: 

GA(t,to)=rexp^*di'P^(i')). 



(7) 



(8) 



Than we do a change of variables P^(t) = G^^{t,to)P^{t)GA{t,to) and D(t) = G^^{t,to)'D{t). 
In this new variables evolution equation M takes form: 



atD(t) = p^(t)D(t). 



(9) 



Now it is easy to see that solution of equation ^ is given by time ordered exponent. We iterate 
this exponent and translate it back to original variables: 



D(t) = GA(t,to)D(to) + 



+ E^=i ( nr=i /to dUeik - h.,)GA{U+i,U)P^{k) ) GA{ti,to)U{to). 



(10) 



For completeness we come back to standard notation which show all the indices explicitly. 
We also perform resummation of virtual part of flavour changing kernel P what leads to 
introduction of flavour changing Sudakov formfactor <^^(t,to) = ft dt'P^^{t'). 



Dk{t,x) = J^ dz'dxoe-'^kit'to)G^^{t,to, z')Dk{to,xo)6{x - z'xo) 



E^=i Efc„_,^...^fc,^A.o nr=i 11 dUeiU - t.-i) /o' dzl !,' dz,e 



-<J>f iti,U-i) 



G'fc,_ifc,_i(*j>*i-l)2i) 



^P^k^.^it^^ ^^)) /o d<+ie-*"(*'*")G4(t, t„, <+i) /o^ dxoDk,{to,xo)5{x - xo DLi z^ m=i 4), 

(11) 
Now to have the full solution we need to find the explicit form of function G^j^{t,tQ, z) which 

represents diagonal transitions. By definition it is a solution of equation ([9]) so we can use 

analogues procedure to solve it. The solution in iterative form is given by: 



G4(t,to,^') = e-*"(*'*«)5(l-z') + 



Eoo 
n=l 



1 ^^/ 



YYl^,JldtMt^-t^-l)Iodz:, 



-*fc(t,tn) 



m=ie-''"^'-''-'^Pg{U,zi) 



s{z' - nr=i zi) 

(12) 
where the initial condition for G^j^{to,to, z') as Dirac delta in z' = 1 was used. 
In this way we obtained hierarchy of two processes. The external- flavour changing process and 
internal - diagonal process. For both processes one can define properly normalised transition 



probability which depends only on previous step (emission) . Probability in the flavour changing 
process is given by: 

(13) 

with normalisation condition: 

/•CXD rl rl 

/ dij V / dzi dz'iUj{ti,ZiZiXi-i,ki\ti^i,Xi^i,ki^i) = 1. (14) 

Jti.i , JO Jo 

For bremsstrahlung process we have: 

Vk{tuZ,z[^^\ti^uz[^^) = e(ti -ii_i)Pfe®fc(ti,Z,)e-*^(*-*-i) (15) 

and normalisation condition: 



/•CO /"l 

/ dti I dZiPk{ti,Ziz[_i\ti-i,z[_i) = l. 
Jti^i Jo 



(16) 



As we can see each of the processes on its own is of Markov type. This means one can use 
standard algorithm to solve each of the two processes separately. What connects this two 
levels are limits for evolution variable. So by means of presented reorganisation'21 we obtained 
hierarchy of two Markov processes which has straightforward interpretation/implementation as 
a generalisation of well known Markovian algorithm. 

3 Results 

Presented algorithm has been used to solve LO DGLAP evolution equation with non-running as- 
It takes into account three massless quarks and gluons. PDFs were calculated for scales reaching 
from Q = IGeV (initial conditions) up to scale of Q = lOOOGeV. Results have been compared 
with standard Markovian algorithm implemented in program EvolFMC ^ ^. Fig. [I] shows gluon 
momentum distribution xDQ{t^x) and quark momentum distribution xDq{t,x) obtained from 
both algorithms and ratio for coresponding distributions for scales equal IGeV, lOGeV, lOOGeV 
and lOOOGeV. Presented results were calculated for statistics equal 2 * 10^ and as we can see 
agreement of both approaches is on the level of 0.5%. 

4 Summary 

The general formalism for reorganization of evolution equation has been used for derivation 
of hierarchical solution to DGLAP equation. Very good agreement with standard Markovian 
algorithm shows that it can compete with it. Generalisation for NLO and for running as case 
can be done and is only a matter of implementation. Also generalisation for different types of 
evolution equations is possible. Especially case of CCFM evolution seems to be worth considering 
because of possibility of simple changing of coupling constants argument for flavour changing 
and diagonal emissions. 
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Figure 1: The left upper plot presents gluon distributions and right upper plot quark distributions obtained from 

hierarchical algorithm (dashed hne) and standard Markovian algorithm-Ev olFMC (straight line) for scales of 1 Gev, 

lOGeV, lOOGeV, lOOOGeV. The lower plots represent coresponding ratio between this two distributions. 
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